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Abstract.  This  paper  is  concerned  with  thermomechanics  of  slender  rods  by  a 
direct  approach  based  on  the  theory  of  a Cosserat  curve  comprising  a one- 
dimensional curve  and  a pair  of  directors  attached  to  every  point  of  the  curve. 

In  all  previous  developments  of  the  thermo-mechanical  theory  of  rods  by  direct 
approach,  only  one  temperature  field  has  been  admitted.  This  allows  for  the 
characterization  of  temperature  changes  along  some  reference  curve,  such  as  the 
line  of  centroids  of  the  (three-dimensional)  rod-like  body,  but  not  for 
temperature  changes  across  the  rod  cross-section.  A main  purpose  of  the  present 
study  is  to  incorporate  the  latter  effect  into  the  theory;  and,  in  the  context 
of  the  theory  of  a Cosserat  curve,  this  is  achieved  by  a recent  approach  of 
Green  and  Naghdi  [1,2]  to  thermomechanics  which  provides  a natural  way  of 
introducing  more  than  one  temperature  field  at  each  material  point  of  the 
curve.  Apart  from  full  discussion  of  thermomechanics  of  rods  and  thermo- 
dynamical restrictions  arising  from  the  second  law  of  thermodynamics  for  rods, 
attention  is  given  to  a discussion  of  symmetries  (including  material  symmetries) 
of  rods  which  in  a reference  configuration  are  straight.  The  paper  also  con- 
tains a detailed  discussion  of  the  linear  theory  of  straight,  elastic,  ortho- 
tropic rods,  including  the  determination  of  the  relevant  constitutive  coefficients. 


^Mathematical  Institute,  Oxford  0X1  3LB,  U.K. 

^Department  of  Mechanical  Engineering,  University  of  California,  Berkeley, 
California  94720,  U.S.A. 
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1. 


Introduction 


This  paper  is  concerned  with  thermomechanics  of  slender  .rods  by  a direct 

approach  based  on  the  theory  of  Cosserat  (or  directed)  curves . A Cosserat 

curve  considered  here  is  a body  ft  comprising  a one-dimensional  curve  (embedded 

in  a Euclidean  3_sPace)  and  two  directors  (i.e.,  deformable  vectors)  attached 

* 

to  every  point  of  the  curve.  The  development  of  a complete  theory  of  a Cosserat 
or  a directed  curve  with  two  directors  begins  with  a paper  of  Green  and  Laws 
[3]  whose  derivation  is  carried  out  mainly  from  an  appropriate  energy  equation, 
together  with  invariance  requirements  under  superposed  rigid  body  motions.  A 
related  theory  of  a directed  curve  with  three  deformable  directors  at  each 
point  of  the  curve,  developed  in  the  context  of  a purely  mechanical  theory  and 
with  the  use  of  a virtual  work  principle,  is  given  by  Cohen  [4].  A further 
development  of  the  basic  theory  of  a Cosserat  curve  along  with  certain  general 
developments  regarding  the  nonlinear  and  linear  constitutive  equations  for 
elastic  rods  is  contained  in  the  more  recent  work  of  Green  et  al.  [5].  For 
clarity's  sake,  we  may  recall  that  the  material  curve  of  ft  can  be  identified 
with  a particular  reference  curve  (often  taken  to  be  an  interior  curve)  in 
the  three-dimensional  rod-like  body,  e.g.,  the  line  of  centroids  of  the  cross- 
section  of  the  rod  in  some  fixed  reference  configuration;  the  directors  at 
each  point  are  regarded  as  representing  the  material  filaments  across  the 
reference  curve,  i.e.,  in  the  cross-section  of  the  rod. 

Throughout  the  previous  developments  of  the  thermo -mechanical  theory  of 
rods  by  direct  approach,  only  one  temperature  field  has  been  admitted  and  this 
allows  for  the  characterization  of  temperature  changes  along  the  reference 
curve  of  the  rod-like  body.  Some  indication  of  how  temperature  changes  across 
* 

The  body  C is  taken  to  model  some  of  the  properties  of  a three-dimensional 
body  of  rod-like  character.  When  the  directors  are  absent  it  reflects  the 
properties  of  a material  curve  appropriate  for  the  construction  of  string 
theory. 


the  reference  curve  of  the  rod-like  body  could  be  dealt  with  has  been  given  in 
the  paper  of  Green  and  Naghdi  [6]  by  using  three-dimensional  approximations , 
but  no  direct  thermo-mechanical  theory  of  rods  with  more  than  one  temperature 
has  been  discussed  in  the  literature  so  far. 

Although  widespread  use  of  the  Claus ius-Duhem  inequalities  has  been  made 
in  three,  two  and  one-dimensional  continuum  thermodynamics,  these  inequalities 
have  been  subject  to  the  criticism  that  in  some  circumstances  they  do  not 
reflect  adequately  ideas  associated  with  the  Second  Law  of  Thermodynamics. 

Green  and  Naghdi  [1]  have  developed  a new  approach  to  three-dimensional 
continuum  thermomechanics  which  is  independent  of  any  particular  mathematical 
expression  of  the  second  law  and  which  imposes  some  restrictions  on  the 
constitutive  assumptions  leading  to  a reduction  of  a number  of  independent 
response  functions  (or  functionals)  in  the  set  of  constitutive  assumptions. 

In  the  present  paper  the  same  approach  is  used  for  a Cosserat  curve  and  this 
provides  a natural  way  of  introducing  more  than  one  temperature  field?  When 
the  directors  are  absent,  the  theory  reduces  to  that  of  a material  curve  which 
may  be  a material  curve  surrounded  by  another  continuum. 

Specifically,  the  contents  of  the  paper  are  as  follows.  Section  2 contains 
a concise  summary  of  the  various  basic  results  of  the  purely  mechanical  theory 
of  a Cosserat  curve  with  two  directors.  With  reference  to  thermal  properties, 
in  section  3 we  admit  at  each  material  point  of  the  curve  of  R a number  of 
different  one-dimensional  temperatures  and  different  one-dimensional  entropies, 
as  well  as  related  thermal  fields;  and,  in  parallel  with  one-dimensional  con- 
servation laws  for  balances  of  mass  and  momenta,  we  postulate  balances  of 
entropy.  Next,  we  recall  the  balance  of  energy  for  the  Cosserat  curve;  and, 
following  the  recent  approach  of  Green  and  Naghdi  [1],  after  elimination  of 

^For  the  purely  mechanical  theory,  it  is  already  clear  how  to  extend  the  theory 
with  more  than  two  directors. 


the  assigned  fields  — i.e.,  assigned  force,  assigned  director  forces  and 
external  rates  of  supply  of  entropy  --  regard  the  resulting  equation  as  an 
identity  to  be  satisfied  for  all  thermo-mechanical  processes.  In  section  4, 
we  briefly  discuss  thermoelastic  theory  of  a Cosserat  curve  on  the  basis  of 
the  new  procedure  in  thermomechanics  (see  section  3)  and  also  compare  the 
results  with  earlier  developments  (see  Green  and  Laws  [3],  Green,  Laws  and 
Wenner  [5])  involving  only  a single  temperature. 

A new  inequality  representing  the  second  law  of  thermodynamics  for  rods 
based  on  the  present  authors'  earlier  work  (Green  and  Naghdi  [1,2]),  along  with 
restrictions  on  heat  flux  vectors  and  the  specific  internal  energy  are  obtained 
in  sections  5 and  8,  respectively,  while  sections  6 and  7 contain  a discussion 
of  relevant  results  for  rods  obtained  from  the  three-dimensional  theory.  The 
last  two  sections,  namely  sections  9 and  10,  are  devoted  to  a discussion  of 
symmetries  (including  material  symmetries)  for  rods  and  the  linear  thermo- 
elastic theory  of  straight  orthotropic  rods.  The  developments  in  sections  9 
and  10  supplement  our  earlier  results  by  direct  approach  (Green  et  al.  [7]) 
for  thermoelastic  rods  in  the  presence  of  a single  temperature  and  previous 
values  for  constitutive  coefficients  in  the  linear  theory. 
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2 . Summary  of  mechanical  theory. 

In  this  section  we  summarize  the  main  kinematics  and  the  basic  equations 
of  the  mechanical  theory  of  rods  based  on  the  work  of  Green  and  Laws  [3]  in 
the  form  developed  by  Green,  Naghdi  and  Wenner  [5]  • A rod  is  a body  ft 
comprising  a material  curve  with  two  deformable  directors  attached  to  every 
point  of  the  curve.  Let  the  particles  of  the  material  curve  of  R be  identified 
with  a convected  coordinate  Q and  let  the  material  curve  in  the  present  configura- 
tion at  time  t be  referred  to  as  c.  Let  r be  the  position  vector  of  c and 

r+S 

d (a  =1,2)  the  directors  at  r.  A motion  of  the  rod  is  then  defined  by^ 


~Qf 


r = r(£,t)  , d = d (£,t)  , [d  cLa  ] > 0 , 


(2.1) 


where 


a3  = a3(C,t)  = dr/3C 


(2.2) 


is  a vector  tangent  to  the  curve  c and  the  directors  d have  the  property  that 

~a 


they  remain  unaltered  in  magnitude  under  superposed  rigid  body  motions.  The 
velocity  and  director  velocities  are  given  by 


v = r(;,t)  , w = d (C,t) 

r*J  /vQf 


(2.3) 


where  a superposed  dot  stands  for  material  time  derivative  with  respect  to  t 
holding  Q fixed.  In  the  reference  configuration  of  R which  we  take  to  be  the 
initial  configuration,  let  the  material  curve  of  R be  referred  to  by  C and 


denote  the  initial  position  vector  by  R,  the  tangent  vector  to  C by  A and 

r>-* 


the  initial  directors  by  D . Then, 


~or 


^The  positive  sign  in  (2.1)3  is  taken  for  definiteness.  Alternatively,  it 
will  suffice  to  assume  that  (d^dgao]  / 0 with  the  understanding  that  in  any 
given  motion  the  scalar  triple  product  [d^d^a^]  is  either  >0  or<0. 
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R = R(C)  = r(C,0)  , A = A (£)  = 3R/dC  = a (£,0)  , 


(2.4) 


D = d (;)  = d (;,o)  . 

~a  ~a 


We  assume  that  the  kinetic  energy  of  the  rod  per  unit  length  of  c is 
given  by 

T =£p(v  • v + 2y°Pv  • w +y“0w  . w ) , (2.5) 

^ ~ -a 

where  p=  p(£,t)  is  the  mass  per  unit  length  of  c and  the  inertia  coefficients 

y°P=  yP°,  y°^=y^®  are  functions  of  £ and  independent  of  t.  We  define  momenta 

per  unit  length  of  c,  corresponding  to  v and  w , as 

~ ~<K 

ST  , . op  . dT  oar  aB  \ 

-=p(v  + ypw^)  , — = p(y  v + y%)  , (2.6) 

respectively.  In  (2.5),  (2.6)  and  throughout  the  paper,  we  use  the  summation 
convention  for  repeated  Greek  indices  over  the  values  1,2. 

With  reference  to  the  present  configuration  at  time  t,  for  each  part  of  c 
between  £=  a,  £=P  We  postulate  the  equations  of  mass  conservation,  momentum, 
director  momentum  and  moment  of  momentum  as  follows: 

r 

a 


dt 


rp 

I p ds  = 0 , 


at  f p(X.+  y° Vds  = f p ~ds  + [~]l  ’ 

a a 

it  f p(y° V y“VdE  = p ic‘ - a33^a)ds  + [£ » 

a a 

it  (I+y°V  + V (y^v  + y^JJds 

or  K 

= f p(r  x f + d X jfca)ds  + [r  x n + d xp°]^  , 

~ ~ ~a  /v  ~a  ~ o' 

where  in  the  above  integrals  the  limits  are  for  values  of  £ equal  to  or  and  (3, 


(2.7) 

(2.8) 
(2.9) 


(2.10) 


as  . «|3a;  . «M  - »3  • >3 


(2.11) 


5. 


and  where  we  have  used  the  notation 


[f(C>t)]e  = f(0,t)  - f(or,t)  . 
a 

Also,  ir.  (2.8)  to  (2.10)  n = n(£,t)  is  the  contact  force  and  pa=pa(r,t)  are 
the  contact  director  forces]"  each  a three-dimensional  field  in  the  present 
configuration;  f=f(g,t)  is  the  assigned  force  and  la - la(Q,t)  are  the 
assigned  director  forces,  each  a three-dimensional  vector  field  and  per  unit 
mass  of  c;  rra=  na(G,t)  are  the  intrinsic  director  forces  which  make  no  con- 
tribution to  the  supply  of  momentum  and  to  the  moment  of  momentum. 

Under  suitable  smoothness  assumptions  the  field  equations  corresponding 
to  (2.7)  to  (2.10)  are 


1 i 


Pa33  " *(0  = P0^33  » 

(2.12) 

an/ac  + \f  = \(v + y0^)  » 

(2.13) 

ap%;  + x£a  = TTa+x(ywWy°^)  , 

(2.14) 

X n + \d^  X (i®  - y°“v  - y°^Wp)  + am/d£  > 

(2.15) 

where 


m = d x p , 


(2.16) 


P = P* ( C ) is  the  reference  density  and  A_  =A_  • A_  is  the  dual  of  (2.11)„. 
oo  j j ~o  ~o  2 


have  the  same  dimensions 


as  r. 


then  n and 


'If 

On  tfre  other  hand,  if  the  directors  are  chosen  to  be  dimensionless,  then  £l 
are  usually  called  contact  director  couples. 


have  the  same  dimensions. 

a 


3-  Thermal  properties.  Thermodynamical  theory  of  rods. 


In  existing  works  on  the  theory  of  a Cosserat  (or  directed)  curve,  only 
one  temperature  field  is  admitted  and  this  is  regarded  as  representing  the 
temperature- variation  along  some  reference  curve,  such  as  the  line  of  centroids 
of  the  cross-section  of  the  rod-like  body.  Also,  the  effect  of  the  thermal 
boundary  conditions  on  the  major  surface  of  the  rod-like  bcdy+  is  incorporated 
into  the  theory  through  the  external  curve  rate  of  supply  of  heat.  The  varia- 
tions of  the  temperature  across  the  rod  cross-section  have  not  been  modelled  so 
far  by  a direct  approach  (within  the  scope  of  the  theory  of  directed  curves ) , 
although  some  indications  of  how  this  could  be  effected  are  implicit  in  some 
work  on  thermoelastic  rods  from  the  three-dimensional  equations  by  the  present 
authors  [6].  As  already  noted  in  §1,  because  of  the  new  approach  to  thermo- 
mechanics of  continua  introduced  recently  by  Green  and  Naghdi  [1],  it  is  now 
possible  to  account  in  a more  general  manner  for  the  thermal  properties  of  a 
rod-like  body  in  a direct  formulation  based  on  a Cosserat  curve. 

Thus,  at  each  material  point  of  the  material  curve  of  R,  we  introduce 

scalar  fields  0=0(£,t)  and  0 =0  (C,t)  (N  = 1,...,K)  representing 

®1®2 " * " ®N  “l0^'"0^ 

the  effects  of  temperature  variation  in  a rod-like  body.  The  curve  temperature  0, 
which  we  require  to  be  positive  (©>0),  represents  the  absolute  temperature  in  the 
curve  c of  the  rod-like  body,  while  the  scalars  © account  for  temperature 

variations  across  the  cross-section  of  the  rod:  the  scalars  0 are 

®l^-'-°N 

assumed  to  be  completely  symmetric  in  the  indices  . . . ,a^.  which  take  the 

values  1,2  only.  Along  with  the  temperatures  © and  0 , we  admit  the 

$ al°i2  “in- 
existence of  external  rates  of  supply  of  heat  r = r(C,t),  r =r  (C,t) 

®1®2 ' ‘ ' °N  °1®2 ’ ' ‘ 

^The  terminology  of  major  surface  refers  to  the  surface  specified  by  (6.4)  in 

§6. 

*The  external  rates  of  supply  of  heat  r and  r ^ include  contributions 

“l“2  • ' • “n 

corresponding  to  heat  fluxes  on  the  major  surfaces  of  the  rod.  They  are  not  the 
same  as  quantities  defined  with  a similar  notation  in  Green  and  Naghdi  [6]. 

7. 


• y-  A A- 


W '* 


1 


per  unit  mass  of  c and  external  rates  of  heat  fluxes  h,  h and 

= = al‘VaN 

h,  h over  each  end  section  of  the  rod.  Also,  we  assume  the  existence 

ala2 ‘ ' ' aN 

of  internal  heat  fluxes  h = h(C»t),  h =h  (r,t)  along  the  rod  at 

ala2---“N  ala2 ‘ * °N 

each  point  £,  in  the  direction  of  increasing  £,  per  unit  length  per  unit  time. 

Each  function  r , . . . ,h  is  completely  symmetric  in  the  indices 

“l<*2 ' • ' °r;  ala2*  ‘ ,aN 

al,a2  ’ ' * ' ,aN*  The  external  rate  of  supply  of  heat  per  unit  mass  of  c is 

defined  as 


r + E r 

a,N  “la2"  ,aN 


(3-1) 


where  the  summation  in  (3.1.  is  over  all  values  of  = 1j2  and  for  all 

N = 1,2,...,K.  Similarly  the  total  internal  heat  flux  at  the  point  Q is  defined 

by 


h + £ h 

a,N  al“2‘ ' ,aN 


(3-2) 

to  tempera - 


We  now  define  the  ratios  of  the  heat  supplies  r and  r 

al°^' ‘ ‘“n 

tures  0 and  0 , respectively,  as  s = s(C,t)  and  s =s  (r, t) 

aiV°fo  aiV“N  aia2'  * ‘°n 

and  call  these  the  external  rates  of  supply  of  entropy  per  unit  mass  of  c. 


Further,  we  define  the  ratios  of  h,h  to  0 and  h , h to 

9 , respectively,  as  the  external  entropy  fluxes  k.k.k 

_ala2 ‘ ‘ ’ ®N  °1VV 


ala2' ' 


k over  the  ends  of  the  curve  c.  Similarly,  we  define  the  ratios  of 

ala2 ’ * ' °N 

h and  h _ tc  the  temperatures  0 and  0 , respectively,  as  the 

r 

= k (C,t)  per  unit 

l2 ' * • 

length  of  c,  in  the  direction  of  increasing  Q.  The  above  definitions  may 
conveniently  be  summarized  by 


to  the  temperatures  0 and  0 

ala2  ,aN 

internal  entropy  fluxes  k = k(£,t)  and  k 


ala2 " ’ * ®N 


r 

i. 


i 


®1«2- 

..Q-;  a1Q^.. 

•“n 

= h 

ala2- 

„ /q 

• ■a.-  a,a0.. 

12 

'°N 

V2--*°N 

= h 

alV 

/e 

• ■ a,-  a-,  a0  • • 

j.  2 

•“N 

ala2 ' ’ ’ aN  * 

k 

“iV 

. k 

• ai;  al“2-- 

V 

finite  limits  as 

9 

0 for 

k = h/9  , k 
k = h/ 9 , k 
We  require  that  the  fields  s 


In  addition  to  the  thermal  fields  already  introduced,  at  each  point  of  the 

material  curve  of  R in  the  present  configuration,  we  assume  the  existence  of 

scalar  fields  ,n  = 1l(C»t)  an(i  H = 11  (C,t)  called  specific  entropies 

ala2  ’ * ’ ®N  “l®2  " ' °t; 

and  internal  rates  of  production  of  entropies  5 = -(£,t),  | = 

“l°2  * ’ ‘ “N 

| (Cj't)  Per  unit  mass  of  c,  where  T1  and  f are  completely 

“lV^N  “iV®!*  «ia2'-'aN 

symmetric  in  the  indices  . . . ,ajj.  The  contributions  of  these  internal  rates 

of  production  of  entropies  to  the  rate  of  production  of  heat  is 

9?  + £ 9 ? (3-4 

^k^-TO-'^n 

per  unit  mass . 

We  postulate  balances  of  entropies  for  every  material  curve  of  R occupying 

* 

a part  j3  in  the  present  configuration  and  write 

« J P^ds  = J p(s+?)ds-  (k]J  , (3-5 

a a 

g 0 

4z  f P T1  ds  = [ p(E  +E  )ds  - [k  . (3.6 

dtJa  1«1VaN  «iV°h  ala2-‘-“N  ala2**'aN0( 

Under  suitable  smoothness  assumptions  it  follows  from  (3*5) > (3*6),  (2.11^  and 

(2.12)  that 


A motivation  for  postulating  (3*5)  to  (3.6)  for  balances  of  entropies  is 
provided  by  consideration  of  derivations  from  three-dimensional  equations 


XT)  = X(b+5)  - Bk/a£  , 


(3-7) 


\T\, 


ala2' ' *®N 


= X(s 


“l«2- 


• ) - ak  /ac 


We  now  introduce  the  first  law  of  thermodynamics  or  the  balance  of  energy 
for  the  Cosserat  curve  ft.  This  states  that  heat  and  mechanical  energy  are 
equivalent  and  that  together  they  are  conserved  for  every  part  of  the  material 
curve  of  ft.  Thus,  with  reference  to  the  present  configuration,  the  balance  of 
energy  may  be  stated  in  the  form 


XT  [ p[e  + i(v  • v + 2y°av  • w + y°®w  . w )]ds 

CLU  J ^ ~ ~ ~ 

= I p[r  + E r + f • v + la  • w ]ds 

a,N  ala2’-,aN  ~ ~ ~ ^ 


+ [n  • v + pa  • w -h-  Eh 

«,N  “lV^V" 


]?.  , 


(3.8) 


where  e = e(£,t)  is  the  internal  energy  per  unit  mass  of  c and  repeated  indices 
are  summed  over  the  values  1,2.  With  the  help  of  (3-7),  (2.12)  to  (2.16)  and 
under  suitable  smoothness  assumptions,  the  field  equation  resulting  from  (3.8) 
is 


X(-e  + 91]  + E 9 71  ) 

a,N  al°(2 ' ’ ‘ °N  al®2'’,0<N 

- X(0?  + E 0 5 ) 

or,N  ^•••ValQ2---°N 

+ p-  k ae/ac  ■ a ae  /ac  = o , 

cr.N  “lV"^  “lV'V 

where  the  mechanical  power  P per  unit  mass  of  c is  defined  by 


(3-9) 


p = aX/a£  + Il*  * ^ + £°  ‘ • (3-10) 

Introducing  the  Helmholtz  free  energy  function  \jr  = \|r(£,t)  per  unit  mass  of  c by 


* = €-971-  E0  71  , 

o,N  ala2“-*N  VV0^ 


(3.11) 


10. 


^ >'V  >-  .•<  ~ 


the  energy  equation  (3*9)  may  be  written  in  the  alternative  form 


+ -ne  + E T\  9 ) 

a.N  ala2‘ ’ ,aN  ala2' ' '“n 

- \(e§  + e e e ) 

cr.N  “lV0^  “lV^ 

+ p-k  ae/dc - Ek  ae  /a;  = o . (3.12) 

a,N  “lQ2,',aN  ala2 ' ' ■ “n 

For  a given  Cosserat  curve  having  a reference  density  p (£),  the  field 

equations  obtained  from  the  integral  form  of  the  conservation  laws  involve  a 

set  of  (K+l) (K+2) + 12  functions.  These  consist  of  the  deformation  functions 

r,d  and  the  temperatures  0,Q  , i.e., 

~a  ‘ " ®N 


[r,d  ,6,6  } , 

~ “lVV 

the  various  mechanical  and  thermal  fields,  namely* 


(3.13) 


,k,k(Y  £V  a.  jv  ,v  } 

~ ~ ~ a±a2  ’ ’ ‘ °N  al°^2  * ‘ ’ “n  al°2 ' ‘ ' °N 


(3.14) 


and 

(f,jfc“,s,s  } . 

We  assume  that  the  fields  (3.14)  are  specified  by  constitutive  equations  which 
may  depend  on  the  variables  (3*13) , their  space  and  time  derivatives,  as  well 
as  the  whole  history  of  deformation  and  temperature.  We  then  adopt  the  fol- 
lowing procedure  in  utilizing  the  conservation  laws*: 

(1)  The  field  equations  are  assumed  to  hold  for  arbitrary  choice  of  the 
functions  (3.13)  including,  if  required,  an  arbitrary  choice  of  the  space  and 
time  derivatives  of  these  functions; 


(3-15) 


*The  density  p is  rot  included  in  (3*14)  and  (3.15)  since,  given  (3.13) , P can 
be  calculated  from  (2.12). 

*For  a more  elaborate  parallel  discussion  in  the  context  of  the  three-dimensional 
theory,  see  Green  and  Naghdi  [1,  §2]. 
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(2)  The  fields  (3- 1*0  are  calculated  from  their  respective  constitutive 


i 


- 

n 


equations ; 

(3)  The  values  of  the  variables  (3-15)  can  then  be  found  from  the  balances 
of  momenta  (2.13)  and  (2.14)  and  balances  of  entropy  (3.7); 

(4)  The  equation  (2.15)  resulting  from  the  balance  of  moment  of  momentum, 
and  the  equation  (3*12)  resulting  from  the  energy  equation,  will  be  regarded  as 
identities  for  every  choice  of  (3.13).  This  will  place  restrictions  on 
constitutive  equations . 

We  note  that  the  quantities  5,5  ,T\,T1  ,i|i  may  be  arbitrary 

“1^2’"%  al°2 ' ’ ’ 

to  the  extent  cf  additive  functions  f,f  ,f,f  , 

ala2'  ’ *®N  aiV'cN 

-0f-  £ 0 f , respectively,  where  f,f  are  arbitrary 

«,N  al<VaN  “l«2---°il  ^ “lV°N 

functions  of  the  variables  (3.13) > their  space  and  time  derivatives  and 
functionals  of  their  histories.  The  additive  functions  have  the  property 
that  they  make  no  contribution  to  the  differential  equations  for  r,d  ,0 
0 and  the  boundary  and  initial  conditions . They  also  make  no  contribu- 

tion  to  the  energy  identity  (3*12)  and  no  contribution  to  the  internal  energy  e. 


We  remove  this  arbitrariness  by  setting 


A.  ^ 

f = f(c) 


A 
= f 


f = 0 


ala2 ' ' * ®N  “la2-’-“N 

, f =0 


(0  . 


(3. 


Then,  the  functions  are  determined  uniquely  and  T|,T1  are 

ala2  ’ " °N  “l'V”0^ 

only  arbitrary  to  the  extent  of  additive  functions  of  Q,  independent  of  t.  The 

in  (3.16)  can  then  be  determined  by  specifying  values 

“l“2*  * *“n 

for  T\,T\ 


A A 

functions  f,f 


*]_«2 ' * ‘ ®N 

in  some  reference  state. 

“l0^  ‘ " °N 

So  far  nc  mention  has  been  made  of  restrictions  on  constitutive  equations 


which  may  arise  from  some  form  of  second  law  of  thermodynamics,  usually 


16) 
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interpreted  in  terms  of  an  "entropy  inequality."  Before  considering  this  and 
in  order  to  gain  some  insight  into  the  nature  of  our  procedure  described  above 


we  study  in  the  next  section  the  relatively  simple  case  of  an  elastic  rod. 

For  later  use,  we  record  the  expressions  for  the  external  work  and  the 
external  heat  supplied  to  any  part  of  the  curve  of  c during  the  time 

interval  t-^StSt  . First,  however,  guided  by  the  results  of  §4  we  observe 

that  in  the  case  of  an  elastic  rod  the  response  functions  ^ ,T|,7\  ,e 

ala2  °N 

depend  only  on  the  vectors  a_,d  ,Bd  /bC  and  the  temperatures  0,0  and 

~3  ~c/  a^.  • .0^ 

are  independent  of  their  rates  and  the  temperature  gradients  d9/d£, 


d9  /d£-  Such  an  elastic  rod  will  be  regarded  as  nondissipative  in  a 

ala2  * ' ’ ®N 

sense  that  will  be  made  precise  later;  and  in  conjunction  with  an  expression 

for  the  external  mechanical  work  supplied  to  any  part  cr  SJS  g of  the  curve  c, 

will  be  used  as  a basis  for  establishing  in  §5  an  inequality  representing  the 

second  law  of  thermodynamics  for  dissipative  materials.  Keeping  this  background 

in  mind,  we  assume  that  the  constitutive  response  functions  for  e,T|  include 

also  dependence  on  the  set  of  variables  a.,d  ,dd  /dC,0,9  ,d9/&C> 

~3  ~0f  aia2 " * ’ °N 

d0  /dC  and  their  higher  space  and  time  derivatives  and  refer  to  this 

“l0^  ’ * ' °N 

set  collectively  as  If.  Further,  let  e'jTl'  denote  the  values  of  e,Tl,  respectively, 
when  the  set  V is  put  equal  to  zero  in  the  response  functions.  Thus,  for  example, 


* H 

li 

1 I 

v 


€ = €(2.3’Sa’aSQ/^’e’eaia2...o^,U)  ’ 


e = e 


* «-J 


(3.17) 


~3  ~a  ®la2 " * ' °N 

v = ,ae/a;,ae  /ac>--«)  > 

~3  ''■or  ~Cr  ala2"""°N 

where  the  dots  in  (3.17)2  refer  to  the  higher  space  and  time  derivatives  of 

a_,...,9  . Then,  with  the  help  of  (2.12)  to  (2.14)  and  the  integral 

al®2  * * ' ®N 
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of  (3.8)  with  respect  to  time,  we  obtain 


ID  = External  mechanical  work  supplied  to  a pan; 
of  the  rod  during  the  time  interval  t^  S 1 5 t 


O $ 

= [ 2|  !*  p(f  • V + la  • w )ds  + [n  • v + pa  • w > dt 


= AK  + AE+Lb+li30 


(3.18) 


and 


W = External  heat  supplied  to  a part  a 5 C-  & °f  the  rod 


during  the  time  interval  t^  5 1 5 tg 

,t2rrP 


= | 2lf  p(r  + £ r )ds  - [h  + £ h ]^V 

Jt1lJa  «,N  0,1q2-“ON  a,N  ala2'  * ‘°N  a 


dt 


(3.19) 


where 


i«t0 


ftp  "B  i*P 

= j j p w ds  dt  , E = J p s ds  , 


dt 


ti  a 


K = ' 


\ p(v  • v + 2y°av  • w + y°^w  • wjds  . 


(3.20) 


The  prefix  & in  (3.18)  denotes  the  difference  operations  on  functions  and  fields 
during  the  time  interval  [t^jt^lj  e.g.,  AE=  E(t^)  - E(t^).  Also,  w in  (3*20)g 


l4. 


is  given  by 


to  = P-  k(e-e')  - \(*,  + T|,e+  Z V 


) 


«,N  ala2 ‘ ‘ aN  ala2 ‘ ‘ ' aN 

- \[('t]-T)/)0  + z (t|  - V )e  ] 

«,N  ala2">aN  ala2*-'aN 

+ \(9|+  Z 9 C ) 

a,N  aia2---ValQ2‘--Q^I 

+ k ae/a;+  z k Be  /ac  , 

a,N  al°i2 ‘ * *N  ala2--‘aN 


t ' = e ' - 0T1 ' - Z 0 


T)  ' 


<*,N  ala2-,,aK 


(3-21/ 


(3-22) 


4.  Thermoelastic  rods 

A thermoelastic  theory  of  rods  by  a direct  approach  was  given  by  Green  and 
Laws  [3]  and  was  developed  further  by  Green,  Lav/s  and  Naghdi  [7]  and  by 
Green,  Naghdi  and  Wenner  [5l-  The  previous  work  made  use  of  a one -dimensional 
Claus ius-Duhem  inequality  and  only  one  temperature  field  was  considered,  which 
corresponds  to  the  curve  temperature  9 of  the  present  paper.  Vie  consider  now 
constitutive  equations  for  a thermoelastic  rod  which  admits  |-(K+l)(K+2) 
temperature  fields  and  we  examine  the  restrictions  imposed  on  these  equations 
by  the  procedure  described  at  the  end  of  §3. 

We  assume  that  the  set  of  variables  (3-14)  are  functions  of 

a,  , d , ad  /ac  > 0 > e , ae/ac  > ae  „ /a;  > (^.i) 

~3  '-a  ~of  c*1a^‘..aN  * aia2‘  ’ ‘°N 

as  well  as  the  reference  values 


h ’ io  ’ ’ 9 


(4.2) 


and  in  addition  may  depend  also  on  the  particle  Q.  In  the  set  of  reference 

values  (4.2),  ® is  the  constant  reference  value  of  9 and  we  have  assumed  that 

the  reference  values  of  9 (a=l,2,;  N=1,2,...K)  are  zero.  Postponing 

ala2  * " ' aN 

the  restrictions  to  be  imposed  by  the  invariance  requirements  under  superposed 
rigid  body  motions  and  recalling  the  procedure  outlined  in  §3,  the  energy  equa- 
tion (3.12)  is  identically  satisfied  for  all  thermomechanical  processes  provided 


a(ae/ac) 


d(ae 7a C 

“iV-ofr 


= 0 , 


n = X > rr“  = \ 

~ 3a,  ~ od 

~3  ~a 


a d* 

1 ' 


(4.3) 


•n  - - & Ti  = . Si 

’ \v“n  aV2...^  ’ 


(4.6) 
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x(9§  + £ 0 % ) + kde/ac+  2 k (ae  /ac)  = o , (4.7) 

«,N  “l«2-”“N  “l“2"*“N  or,N  al°^  * ’ " aN  “iq2'*-“n 

where  as  indicated  in  (4.3),  the  function  i|j  is  independent  of  temperature 

gradients.  It  is  understood  that  Q appears  in  * only  in  a form  which 

• • • Qfjj 

is  completely  symmetric  in  the  indices  Formulae  (4.5)  end  (4.6), 

with  9 absent,  were  obtained  by  Green,  Naghdi  and  Wenner  [5]"'"  with 

“l0^  • • • Ofo 

the  help  of  the  Claus ius-Duhem  inequality.  Alternative  forms  for  the  results 
(4.4),  (4.5)  and  (4.6)^  have  been  given  by  Green  and  Laws  [3]  and  by  Green, 

Naghdi  and  Wenner  [5]  after  making  use  of  invariance  conditions  under  super- 
posed rigid  body  motions.  For  later  use  we  record  one  of  these  alternative 
forms  here  and  for  this  purpose  we  introduce  the  notations 


d_  = a , h.  . = d.  • d.  , \ . = d.  • ad  /dC 
~3  ~3  ij  ~i  Qfi  ~i  w 


i1  • d1  = 61  , hlJ  = d1  • dj  , h = 


det  h. 


ij  ’ 


where  Latin  indices  take  the  values  1,2,3  and,  from  (2.1),  we  recall  the 
condition 

hi  ■ ISaA.^)  > 0 * C*- 

1 

so  that  h2  is  a single-valued  function  of  h.  ..  For  the  duals  of  quantities  in 
(4.8)  defined  relative  to  a reference  configuration,  we  employ  the  corresponding 
capital  letters,  e.g.,  A^,  H^j,  etc. 

Under  superposed  rigid  body  motions  the  vectors  a3’da’ddc/®C  tecome 
Q ~3’  §,  where  Q is  a proper  orthogonal  tensor  function  of  the 

time,  and  the  corresponding  value  y+  of  the  free  energy  response  is  given  by 


In  this  paper  p.  490,  equation  (7*5)^  should  read  £ = dd 


**  - *(«l3,9ao.aa4/ac’9’9«1ofe...^!Ss>S.>iS0/!l£’eiC) 

Then,  by  (4.4),  (4.10)  and.  the  invariance  condition  v = we  have 


(4.10) 


(4.11) 


for  all  proper  orthogonal  tensors  Q.  It  follows  from  Cauchy's  representation 
theorem  that  % may  be  expressed  as  a different  function  of  the  inner  products 
and  scaler  triple  products  of  namely  the  inner  products  (4.8)^  ^ and 


(4.12; 


[d.ad^/ac  a^/ac!  . 


In  view  of  (4.9),  each  of  the  scalar  triple  products  in  (4.12)  may  be  expressed 

l 

as  a single-valued  function  of  h2  and  (4.8)2  3 En^  hence  °f  and  ^ai' 

Similarly,  if  instead  of  (2.1)-  or  (4.9)  we  make  the  choice  [d.  ,do,d_]<0  for 

3 ~l  ~3 

all  motions,  tnen  we  may  again  reduce  to  depend  only  on  h^A^,  aPai"t  from 
the  reference  variables  and  temperatures.  Hence,  we  may  replace  (4.4)  by  the 
different  single-valued  function 


(4.13) 


where 


vij  hij  ‘ Hij  ’ Kai  \*i  ‘ Aoi  * 


(4.14 


Let  a.  be  a set  of  base  vectors  with  fa. a„a_]j(0  and  let  aE  be  their  reciprocals. 
The  base  vectors  a^  ordinarily  will  be  taken  to  be  orthogonal  but  this  is  not 
essential  at  this  point  in  our  development.  It  is  now  convenient  to  introduce 
the  component  forms  of  the  kinematic  variables  cL  ,d^  and  of  the  kinetic  variables 
relative  to  the  base  vectors  a.  or  a^.  Thus,  we  write 


* s X ' 


(4.15) 


r 


d.  = d.,aJ  = dA,  , d*  = d^a,  = dA'3  , 

~i  ij~  i-~j  ~ ~J  j~  ’ 


n = n a. 

~ ~i 


a oa  a al 

TT  = TT  8 . , P = P a 


Si  • 


It  then  follows  from  (4.5),  (4.6)  and  (4.13)  that 


d3  v . k v1)  - sx  -*■*- 


dy 


33 


, d»y-d3;x  “pvl.x-jt-  , 

i i v 3Y 


d8;(,o,i-xvy1)*d«:(„»1-ivy1)  - 4X  j*-  , 


o9 


(4.16) 


g- 

o a«oJ 


A 

M. 
‘ de 


A 

.14- 


ala2 ’ ‘ * ®N 


ae 


ala2---aN 


nd 


and  in  evaluating  (4.16),  _ , , ^ is  regarded  as  a function  of  -y  ,y  _ a 

-*->^>3  33  cr3 

?(yap  + ,Ypa)-  The  results  (4.16)  are  equivalent  to  (7-40)  in  Green,  Naghdi  and 
Wenner  [5];  those  in  Green  and  Laws  [3]  are  a special  case  of  (4.16)  when 
d.  = a.  , apart  from  some  changes  in  notation  and  definitions. 

With  the  help  of  (2.14)  and  (2.16),  the  moment  of  momentum  equation  (2.15) 
can  be  expressed  in  the  form 


&3*£.+Ve“+  (**J*u*za  = 0 • 


(4.17) 


This  equation  is  identically  satisfied  by  the  expressions  (4.16).  This  implies 
that  the  only  relevant  field  equations  are  (2.13)  and  (2.14). 

In  a similar  manner,  with  the  help  of  invariance  conditions  under  super- 
posed rigid  body  motions,  the  entropy  flux  functions  k,k  may  be  reduced 

al°2‘  ”®N 

to  depend  on  the  variables  displayed  in  (4.13). 

With  the  help  of  (4.5)  and  (4.6)  we  see  that  the  expression  for  w in 
(3.2l)  is  zero  and  the  external  mechanical  work  Ifl  and  external  heat  supplied 
ft  in  (3.l8)  and  (3.19)  reduce  to 


19. 


;r 


16=  AK+  AE'  + ID 


(4.18) 


I 


5.  The  second  law  of  thermodynamics 

Previously,  in  the  context  of  the  three-dimensional  theory,  Green  ar.d 
Haghdi  [1,2]  have  discussed  the  nature  of  thermodynamic  irreversibility 
arising  from  a mathematical  interpretation  of  a statement  of  the  second  lav 
namely  that  ”it  is  impossible  completely  to  reverse  a process  in  which  energy 
is  transformed  into  heat  by  friction."  Here  we  follow  the  same  procedure  and 
reconsider  a mathematical  interpretation  of  a second  law  appropriate  for  a 
direct  theory  of  rods  which  admits  more  than  one  temperature  field.  In  earlier 
works  [3,5]  > a Claus ius-Duhem  inequality  was  used  and  only  one  temperature 
field  was  admitted. 

The  state  of  a rod  at  time  t which  is  regarded  as  representing  a thin 

rod-like  three-dimensional  body,  is  described  by  the  position  vector  r and  the 

directors  d , the  velocities  v,w  , the  temperatures  9,0  (N = 1,2, . . . ,K) 

~o  ~ al°2  ‘ 

throughout  the  curve  c,  together  with  the  constitutive  response  functions  for 

the  fields (3. l4) . A thermo-mechanical  process  or  simply  a process  is  a time 

sequence  of  states:  it  is  a continuous  oriented  curve  in  the  space  of  states, 

i.e.,  the  (9,9  ,h..,X.  ) space.  Thus,  a process  may  be  specified  by  a 

al°^ ' * • ®n  ^ 0(1 


sequence  of  values  of 


r , d ,9,9 
~ ~a  q'1®2  ' ' ’ ®N 


(5-1) 


on  c in  the  time  interval  OStSa.  Similarly,  the  reverse  process  is  a process 
defined  by  a sequence  of  values  of  (5.1)  on  c in  the  time  interval  (jStS2c 
subject  to  the  conditions 


r(t)  = r(2o-t)  , d (t)  = d (2a-t)  , 

'X'  ~ r-wQ;  *xQ( 

0(t)  = 0(2o-t)  , e „ (t)  = 0 „ (2a- t)  . 

a1a2‘ ' *®N  “l®2,,,0N 

In  any  process  the  work  done  by  the  external  mechanical  forces  acting  on 


(5.2) 
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a part  a5£S0  of  c,  and  given  by  (3.18),  is  a positive  or  negative  depending 
on  whether  the  external  work  is  supplied  to,  or  is  withdrawn  from,  the  part. 

In  general,  some  of  the  work  done  results  in  a change  of  the  kinetic  and  internal 


energies  represented  by  the  first  two  terms  on  the  right-hand  side  of  (3.18)  , 

each  of  which  may  be  positive,  negative  or  zero.  Also  part  of  the  work  done 

may  be  positive  with  a corresponding  absorption  of  heat  by  the  part 

or  negative  with  a corresponding  absorption  of  heat  by  the  part.  We  note  that 

in  the  case  of  an  elastic  material  the  different  contributions  to  Id  will  vary 

in  sign  depending  on  the  process  and  will  not  be  restricted  to  be  either 

positive  or  negative  for  all  processes.  Consider  any  process  in  the  time 

interval  OstSa  and  its  reverse  process  in  the  time  interval  <j§tS2 a.  If 

the  process  is  reversed  in  such  a way  that  at  the  end  of  the  process  and  its 

reverse  process  the  elastic  rod  has  returned  to  its  original  state  with  A0=O, 

A0  = 0,  Ah, , = 0,  AX  , = 0,  Av=  0,  Aw  =0  and,  hence,  Ae  = 0,  ATI  = 0, 

al°^  ’ ' ' Ofl  ~ ~ '"■Of  '** 

45«Lv  V°’  ‘k“°« 

Ak  =0  and  AK=0,  AE=0,  then  all  the  work  done  in  the  process  is 

recovered  as  work  in  the  reverse  process.  This  recovery  of  work  would  not  be 
possible  if  in  every  arbitrary  process  part  of  Id  always  has  a positive  sign. 
With  this  motivation  in  mind,  we  assume  that  for  any  arbitrary  process  in  a 
dissipative  rod  only  part  of  the  work  done  is  recoverable  as  work  in  the 
reverse  process,  the  rest  being  transformed  into  heat.  We  therefore  assume 
that  in  every  process  part  of  the  work  done  is  always  nonnegative.  Then,  if  at 
the  end  of  any  process  and  its  reverse  process  the  rod  has  returned  to  the 
same  state,  some  of  the  work  done  is  always  transformed  into  heat.  Recalling 
that  Ui,  = 0 in  (3*18)  in  the  case  of  an  elastic  rod,  we  interpret  the  above 
assumption  for  a dissipative  rod  by  requiring  that 


‘if  work  is  extracted  in  the  process  then  it  is  absorbed  by  the  rod  in  the 
reverse  process. 
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' '***«*r 


(5.3) 


te2  5 o 

for  all  parts  of  c and  all  processes,  where  1&,  is  given  by  (3. 20).  Since  t^,t 
are  arbitrary  and  w has  already  been  assumed  to  be  continuous,  it  follows  that 


\w  = P-  \(e-e  ')  - + T]'0+  E T\'  „ 9 ) S 0 

«,N  VV°N 

for  all  thermo-mechanical  processes.  Also,  from  (3-19)  and  (5.4),  we  have 

pP  • . 

M5  \( 9T1 ' +29  T|'  )ds  dt 

*^t1  a a,N°'l^"'Q!Nala2"'0'N 

so  that  the  external  heat  supplied  to  any  part  of  c is  bounded  above  in  any 


(5-4) 


process. 


l 


6.  Summary  of  results  from  three-dimensional  mechanical  theory. 

Consider  a three-dimensional  body,  embedded  in  Euclidean  3-space,  and  let 

• i / » * 

its  particles  be  identified  by  convected  coordinates  £ (i  = 1,2,3).  Let  r 
denote  the  position  vector,  relative  to  a fixed  origin,  of  a typical  particle 
of  the  three-dimensional  body  in  the  present  configuration  at  time  t.  Then, 

r = r (£  ,£  ,£  ,t)  , ^ = 3r  /a;  , V = r , 

« 

gik  = gi  ' Sk  ’ gl  • & = 6k  ’ glk  = gl  • £ ’ g = det  §ik  ’ 

where  g.  and  g1  are  coveriant  and  contravariant  base  vectors,  respectively, 
g and  g are  covariant  and  contravariant  metric  tensors,  respectively,  and 

IK 

6^  is  the  Kronecker  delta.  Also  a superposed  dot  denotes  material  time  deriva- 
K 

i * 

tive  holding  £ fixed  and  v is  the  velocity  vector. 

The  stress  vector  t across  a surface  in  the  present  configuration  whose 
* 

unit  outward  normal  is  v may  be  expressed  in  the  form 


* i i k * ik 


, -_i  / 2 " J-n.  * * i *i 

t = ViT  /g  = VfT  & , V = vtg  = V g.  , 


(6.2) 


where  t are  the  contravariant  components  of  the  symmetric  stress  tensor. 

We  do  not  recall  here  the  consequences  of  the  conservation  laws  of 

the  three-dimensional  theory  since  they  will  not  be  needed  in  the  present 


paper . 


The  parametric  equations  £ = 0 define  a curve  in  space  at  time  t which  we 

assume  to  be  smooth  and  which  we  identify  with  the  curve  c.  Any  point  of  c is 

specified  by  the  position  vector  r relative  to  the  same  fixed  origin  to  which 
* 

r is  referred,  where 


r = r(£,t)  = r (0,0, £,t)  , £'  = £ . 


(6.3) 


Let  the  boundary  of  the  three-dimensional  region  occupied  by  the  body  at  time  t 
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be  specified  by  the  material  surface 


(6.4) 


and  by  the  surfaces  £ = <2,9.  The  surface  (6.4)  is  such  that  £ = constant  are 
curved  sections  of  the  body  bounded  by  closed  curves  on  this  surface. 

Suppose  now  that  r in  (6.1)  is  a continuous  function  of  £“,t  and  has 
continuous  space  and  time  derivatives  of  order  2 in  the  bounded  region  lying 

inside  the  surface  (6.4)  and  between  £=a,9-  Hence,  to  any  required  degree  of 

* 12 
approximation,  r may  be  represented  as  a polynomial  in  £ ,£  with  coefficients 

which  are  twice  continuously  differentiable  functions  of  £,t.  Instead  of  con- 
sidering a general  representation  of  this  kind,  we  restrict  attention  here  to 
the  approximation 


r = r(£,t) + £ad  (£,t)  . (6.5) 

~ ^ '■'Of 

Given  the  approximation  (6.5)  it  is  known  (see,  e.g.,  Green,  Laws  and  Naghdi 
1968,  Green  and  Naghdi  1970)  that  the  field  equations  of  the  forms  (2.12)  to 
(2.15)  can  be  derived  from  the  three-dimensional  equations  provided  we  identify 
d in  (6.5)  with  (2.l)0  and  adopt  the  definitions 

~Q  2 


pa33  = X = LUdC1(i£2  ’ n = P*S^  = n(C:>C2)  > 

4y“6  = Jj  M,  C^d^2 


(6.6) 


2,  = JJ  I3d£ld£2  > IT*  - JJ  Tad£xd£2  , 
E“  - JJ  £“E3d£V  , 


(6.7) 


where  p is  the  three-dimensional  mass  density  and  the  integrals  are  taken  over 
any  surface  £= constant  bounded  by  (6.4).  For  some  purposes  it  is  convenient 
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In  this  section  we  obtain  seme  thermodynamical  results  for  a rod-like  body 
on  the  basis  of  the  recent  thermodynamical  theory  of  Green  and  Naghdi  [1], 

Thus,  along  with  the  three-dimensional  temperature  field  9=9  ( ,-b ) > 0 we 


admit  the  existence  of  an  external  rate  of  supply  of  heat  -h  per  unit  area 

* * 

acting  across  the  boundary  SR  of  a region  of  space  R occupied  by  the  body  in 

the  present  configuration  at  time  t.  Also  we  assume  the  existence  of  an 

• v * / i * 

internal  surface  flux  of  heat  -h  = -h  (£  ,t  ; v ) per  unit  area  across  each 

, . * * * 
sun  ace  3P  which  is  the  boundary  of  an  arbitrary  part  p of  R . We  define  the 

ratio  of  the  heat  supply  r to  temperature  9*  as  s*=s*(£i,t)  and  call  this  the 

external  rate  of  supply  of  entropy  per  unit  mass.  Similarly,  we  define  the 
— * * 

ratios  of  h and  h to  temperature,  respectively,  as  the  external  rate  of 
surface  supply  of  entropy  k per  unit  area  of  dR  and  the  internal  surface  flux 
of  entropy  k =k  (r^,t  ; v*)  per  unit  area  of  dp  • Thus 

* * * _*  *_*  * * * , 
r = 6 s , h = 9k  ,h=0k  . (7-1) 

-x-  * i 

In  addition,  throughout  R we  assume  the  existence  of  a scalar  field  T)  = T)  (£  ,t) 

per  unit  mass,  called  the  specific  entropy  and  an  internal  rate  of  production 
* *-  i x 

of  entropy  £ = £ (C  >t)  per  unit  mass.  The  contribution  of  the  latter  to  the 

* * 

rate  of  production  of  heat  is  0 5 per  unit  mass. 

We  recall  the  balance  of  entropy  in  the  form 


for  every  material  volume  occupying  a part  p in  the  present  configuration.  It 


follows  that  k is  linear  in  v 


where  p is  the  entropy  flux  vector.  Then,  from  (7.1)  and  (7.3),  h 


and  we  may  define  the  heat  flux  vector  q by 


* * * 
q = 0 p 


(7.4) 


With  the  help  of  (7-3)>  the  field  equation  corresponding  to  (7.2)  is 


*.  * * . * ■*. 


p T)  = p (s  +|  ) - div  p , 


(7.5) 


where 


div  p*  = g‘?d(p*kg2)/d;k 


A “2  A 


(7.6) 


Multiply  (7.5)  by  £ £ ...£  and  integrate  over  an  arbitrary  part  p in 


the  present  configuration.  After  using  (7.3)  and  some  straightforward  manipula- 
tion, we  obtain 

d f * * al  a2  aN  P */  * *1  a2  aN  !*  * ai  °b  “n 

dtj.pTl£  c -C  ( •••£  • #k£V...;V 

p p Jdp 

+ J *(P  1£^...£Q!n+ ••• +P  ^£ai...£Vl)dv  . (7.7) 

JP 

* 

We  now  suppose  that  P is  a region  bounded  by  the  surface  (6.4)  and  by 
£ = o,p.  Then  for  a rod-like  body,  from  (7*2)  and  (7*7)  we  may  derive  the 
balance  equations  (3*5)  and  (3.6)  without  any  approximations  provided  we  make 


the  following  identifications: 

[•  * 1 ? I*  * al  a?  aTJ  l ? 

XT]  = J nT|  d£  d£  , XT)  a = p-l)  £ £ •••£  "dpiC  , 

J “la2 ' * ' aN  j 


..  T *AA  .Y.1.,2 

4s  jv  ~ - J £ ,,  ...  £ d£  d£ 

ala2  aN  J 


(7.8) 


Xs  = jj  tis*d£ld£2-|  k"g2[(v*1-X1vX3)d£2-  (v^-X2v*3)d£1]  , (7 -9) 


- I ®Nr,  *1  ,1  , *2  ,2  *3\J„11  fr, 

j ^ 8 = t ...£  L (v  -X  v )d£  - (v  -X  v )d£  ] , (7.10) 


H = JJ  M-?*d£1d£2  , 


(7.11) 


' I 

i 


9 

A 

4,  I 

| 


*5 


a,  a0  • • - a.. 

J.  c.  *\ 


_ f*  ~ * ai  ®N1P 

XV«V«!i  + JJ  C -C  dC  dC 


(7.12) 


*°foA  .*t:-i>,-i,-2 


r«  1 '-<*1  <*p  a~  a,T 

" JJ  g (p  c C -‘-C  +-+P  )iC*C  > 


J k*(^AqW  - Jj  p*3AcV  , 

. -ir^g^^.-cVK? 

> • * 0 w 

' -JJp*W1...c“NdcV  • 


a-  a. 


(7.13) 


We  also  recall  the  three-dimensional  equation  for  the  conservation  of 
energy,  namely 

dt  J • v + e ) p dv  = J #(s  e + f • v )p  dv  + | *(t  • v - k Q )da  , (7.1U) 

P P dp  ~ ~ 

* 

where  e is  the  internal  energy  density.  Suppose  in  addition  to  the  approxima- 
tion (6.5)  for  the  displacement  vector,  we  adopt  the  approximation 

a.. 


0 = 9+  E 9 „ 

a,N  “lV  • ,aN 


cV2.*.fN  (9>0) 


(7.15) 


for  the  temperature  field.  The  summation  in  (7.15)  is  over  all  values  of 

N = 1,2,...,K  and  all  Greek  indices  have  the  values  1,2.  The  functions 

9,9  depend  on  C,t  and  0 is  completely  symmetric  in  its 

“lV^  “lVaN 

suffices.  Then,  for  a rod-like  region  bounded  by  the  surfaces  (6.4)  and  £=<*,(3, 
from  the  energy  equation  (7*14)  we  can  derive  the  equation  of  balance  of  energy 
(3.8)  for  a directed  curve  provided  we  make  the  identification 


T * 1 2 

X«  = jj  ne  d£  d;  , 


(7.16) 


in  addition  to  (7.8)  to  (7.13). 
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8.  Restrictions  on  heat  flux  vectors  and  internal  energy 

Suppose  that  the  three-dimensional  rod-like  body  is  in  equilibrium  with 

* 

v = o and  all  response  functions  are  independent  of  the  time.  As  part  of  their 
thermodynamic  restrictions  on  constitutive  equations,  Green  and  Naghdi  [1] 
have  adopted  the  classical  inequality 

* *•  * * 

-q  • g 5 0 or  -p  • g 5 0 , 


£ = grad  6 


for  all  time-independent  temperature  fields.  Integrating  the  inequality  (8.1), 

c 

1 2 

with  respect  to  £ and  Q , we  obtain 


ni  * * 1 ? 

g*£  . £ dC  at  « o . 


(8.2) 


With  the  approximation  (7-15)  for  0 and  with  the  help  of  (7.12),,  and  (7-13) , 
it  follows  from  (8.2)  that 

x.  e 0 ? „ +kde/ac  + l k ae  /ac  s o (8.3) 

cr,N  al°fe* * *®N  *1*2’ * *aN  <y,N  “lV^N  “lV^N 

for  all  equilibrium  displacement  and  temperature  fields.  With  the  above  motiva- 
tion, we  add  an  inequality  of  the  form  (8.3)  for  all  equilibrium  states  to  the 
thermodynamic  inequality  (5-4)  which  was  derived  directly  from  one-dimensional 
postulates . 

Next,  suppose  that  the  Cosserat  curve  ft  is  at  rest  with 

v = o , w =0  (8.4) 

for  all  time  and  with  the  deformation  gradient,  director  and  director  gradient 

each  constant  for  all  time.  Then,  by  (2.12),  p is  independent  of  t.  In 

addition,  we  restrict  the  temperature  fields  to  be  spatially  homogeneous  so 

that  0=9(t),  9 =0  (N=1,...,K).  Keeping  these  conditions  in  mind, 

®L°fe"  ,0fN 

from  a combination  of  (3-7)  and  C3 - 9 ) we  have 
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(8.5) 


In  the  rest  of  this  paper  we  restrict  attention  to  the  case  in  which  K = 1 
in  the  earlier  part  of  the  paper  (§§3»4,5)  so  that  we  consider  only  three 
temperature  fields  9,0^, 0£.  The  inclusion  of  the  two  temperatures  0^,0^ 
allows  us  to  take  some  account  of  temperature  variations  across  the  thickness 
of  the  rod  but  the  more  general  situation  K>1  can  be  dealt  with  in  a similar 
way.  We  also  restrict  attention  to  rods  which  in  a reference  configuration 
are  straight.  We  therefore  consider  the  form  of  the  Helmholtz  free  energy 
function  in  (4.13)  which  is  such  that  the  rod  described  in  §2  models  the  main 
features  of  a three-dimensional  elastic  rod  which  has  the  following  properties 
in  its  reference  configuration: 

(i)  It  is  straight  with  constant  cross-sectional  area6  normal  to  the 
line  of  centroids  of  sections,  (ii)  each  cross-section  has  two  orthogonal 
axes  of  symmetry  through  its  center  of  mass,  (iii)  it  is  homogeneous  and  of 
constant  temperature,  (iv)  the  material  of  the  rod  possesses  rhcmbic  symmetry 
with  respect  to  orthogonal  axes  of  symmetry  and  the  line  of  centroids. 

We  choose  the  initial  reference  curve  C to  be  along  the  direction  of  the 
constant  unit  vector  e_,  and  the  initial  directors  in  their  reference  con- 
figurations  to  be  specified  by 
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state  has  a position  vector  R specified  by 

l - tZ-3* ■ 0-2) 

Materials  with  rhombic  symmetry  with  respect  to  the  three  directions  e.  contain 

~i 

three  classes,  namely  rhombic-pyramidal,  rhombic-disphenoidal  and  rhombic- 
dipyramidal.  For  each  of  these  classes,  it  can  be  shown  that  the  three-dimensional 
energy  function  for  an  elastic  material  is  also  form- invariant  under  the  three 
separate  coordinate  transformations 

(a):  C1--;1  , (b):  Q2  -- f , (c):  C— C • (9-4) 

The  conditions  (9-4)  are  those  which  arise  from  symmetry  restrictions  usually 
called  orthotropic  symmetry.  In  the  direct  formulation  of  the  theory  of  rods 
under  discussion,  we  replace  the  notion  of  rhombic  symmetry  by  that  of  orthotropic 
symmetry  with  respect  to  three  orthogonal  vectors  e. . Therefore,  we  further  assume 
that  the  one-dimensional  energy  function  (9.2)  is  unaltered  in  form  under  each  of 
the  transformations  (a)  to  (c)  listed  below: 


(a): 


D1  ‘‘-“I  ’ D2  -D2  ’ *L  — £l  ’ &>  -&»  ’ 

® “*  ® » ^2.  ” ®i  5 02  02  ’ 


i.e.,  under  the  transformation 


^1  ^Y12,Y13’^12,^21,^13,01’^1^  ”*”^1  ’ 


Y11  ^ Y11 

t 

CVl 

& 

y22  ’ 

y33  y33 

9 

y23  “*  y23  ’ 

K11  *11 

’ K22  "* 

K22  ’ 

K23  “■  K23 

9 

9-9  , 92  - 92 

Do  - - D0 

- D , d 

— * - 

d.  « d.  -♦  d • 

2 2 

* 1 

1 ~2 

* it 1 ita.  * 

0>):  j 

9 — 9 , 

82  * ®1 

0l 

i.e.,  under  the  transformation 


+See,  for  example,  Qreen  and  Adkins  [8]  which  contains  additional  references 
on  material  symmetries. 
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^2  lv^2  * ^12 9 ^21 9 ^2*  * ^2 ^2  9 


and 


(c): 


'll  -*  Y11  ’ 

Y22  "* 

y22  ’ y52 

“ v33 

:11  ■*  K11  » 

*22  "* 

*22  ’ *13 

- ic13 

- D.  , D_ 

- D„ 

, D_  — - D_ 

, d. 

1 ’ 2 

2 

’ ~6  ~3 

"*  9 , 9X  - 

91  ’ 

9g  — 92  j 

(9-6) 


9-0  , 9i  "*  ei 


’(9*7) 


i.e.,  under  the  transformations 


y13  v13  ’ Y23  Y23*  *11  <11  ’ *12  *12  * *21  *21  ’ *22  *22  ’] 

Y11  “*  Y11  * y22  " y22  ’ Y33  ■*  Y33  ’ *13  “ *13  ’ *23  ^ *23  ’ 

9 — 9 , 9-^  — 9^  > 9g  “•  02 

Next,  in  view  of  the  condition  (ii),  we  assume  that  in  (9*2)  is 

unaltered  in  form  under  the  transformation 


(d) :  £i  ~«i  ’ ’ ~3  ”*  ^3  ’ 

i.e.,  under  the  transformations 

Dl'*'Dl  ’ D2  " D2  ’ Yij  **  Yij  ’ *icr  *icr  * 9 “*  9 ’ 9cr  " 9a  ’ (9,8a) 

and  is  unaltered  in  form  under  the  transformation 

(e) :  *1  - Si  » *2—22  > % - % » 

i.e.,  under  the  transformations 

D1  “*  D2  » °2  D2  » Yij  "*  Yij  » *ia  “*  *icr  » 9 - 9 » 9o  ^ 9a  ‘ (9>8b) 
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^ ‘ - r*  ■ ' & ^j*f*W*" 


M 

r 


'Vg-rr  -r- 


j 


A similar  discussion  may  be  carried  out  for  the  functions  ,k,k 

a a 

-xoep'  that  each  of  these  functions  may  depend  on  3©/3G>  d8./d£  in  addition  to 

P 

these  specified  in  (9-2).  In  view  of  the  energy  identity  (3-12)  we  assume 

that  *9,  e 0 , k39/3C,  k 39/dC  for  a =1,2  are  each  unaltered  in  form  under  the 

3 a a 

transformations  (9.5)  to  (9.8).  To  each  of  the  transformations  (9-5)  to  (9.8) 
ve  add,  respectively,  the  transformations 

(a) :  a9/ac  - ae/ac  , aej/ac  -- a^/ac  > ae2/ac  - ae2/ac  , 

(b) :  a©/a;  - ao/ac  , a0x/ac  - a0x/ac  , ae2/a;  — ae2/a;  , 

(9-9) 

(c) :  ae/a; - ao/ac  , ae^ac -- ao^/ac  , 

(d)  and  (e):  39/3;  - ae/ac  , a© ^ac  - a© yac  . 

After  applying  invariance  conditions  to  ^,§,|a>k,k^  using  (9*5)  to  (9*9)  we 
specify  the  initial  directors  to  be  coincident  with  the  unit  vectors  e,,e„, 
i.e. , 


= 1 > d2  = 1 . (9.10) 

In  order  to  make  the  above  conditions  explicit  we  limit  our  attention  to 
the  linear  theory  of  a thermoelastic  rod  which  is  in  equilibrium,  unstressed 
and  at  uniform  temperature  in  its  reference  configuration.  For  this  linear 
theory  the  position  vector  r,  directors  d and  their  corresponding  velocity 
fields  assume  the  forms 


r = R + u , d.  = D . + 6. 

f+*  rsj±,  r*J±.  i 


u 


ui£l 


Ui*i 


w 

~or 


5 .e. 


(9-11) 


and  we  replace  9 by  ®+0  so  that  are  small  and  © is  the  constant 

reference  temperature  of  the  rod.  Using  (9. 10),  the  kinematic  measures  of 


deformation  become 


I 


v.  . = e . • 6 . + e . • 6.  - 6.  .+  6..  , 

Yij  ~i  ~j  ~i  ij  Ji 


ij  JJ 


(9-12) 


Kai  = a6cri/3£  » 63i  = • 


The  response  functions  n,rr  ,p  and  external  fields  are  referred  to  the 

orthonormal  basis  e. . Thus 


n = n . e . , tt  = rr  . e . 

1~0_  ''■Qf  <2l~l 


’ £a 


P -e- 
ai~a. 


f = f.e.  , l = l .e, 
.w  i~i  -*qj  ai~i 


(9.13) 


and  the  relations  (4.16''  simplify  to 


A A 

n = 2p  , n = p r^- 

3 3y33  « dVa3 


A 

ti  = 4p  dt 

aP  Pa  aYag 

"‘-le  > \ 


. D =Q  ' 

’ ^i 


A 

. *L 

d9a 


(9.14) 


where  p is  now  the  constant  reference  density  and  in  evaluating  (9-l4)0  _ the 
A ^ ><3 

response  function  \|r  is  regarded  as  a function  of  and  + The 

equations  of  motion  (2.12)  to  (2.15)  are  replaced  by  their  corresponding 

linearized  forms  in  terms  of  n.  ,tt  . and  p . with  X=  p,  and  the  linearized 

i ori  oi  K 

entropy  balance  equations  still  have  the  forms  (3*7)  with  \=  p.  All  these 
equations  are  listed  below  in  four  groups. 

We  omit  detailed  discussion  of  the  invariance  conditions  satisfied  by 
which  is  now  a quadratic  form  in  note  that  the  final  results 

separate  into  four  groups  corresponding  to  flexure  (2  modes),  extension  and 
torsion  of  the  rod.  Similarly  the  symmetry  restrictions  on  §,§a,k,k^,  which 
are  now  linear  forms  of  degree  one  in  Y^,Kai,9)9a>d0/dC  and  30 place  these 
functions  in  the  same  four  categories.  In  addition,  these  functions  are  subject 


i 

> 

i 
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to  the  identity  (4.7) > which  in  the  present  context  becomes 


p (0+0)5  + pe151  + pe2i2  + kae/3;  + + k2ae2/d£  = o , 


and  the  inequality  (8.3) » namely 


+ p02?2  + kae/dc  + k1ae1/sc  + k2&e2/dc  s o 


(9.16) 


Symmetry  restrictions  of  the  type  considered  in  this  section  must  also  apply 
to  the  kinetic  energy  (2.5)  so  that  the  inertia  coefficients  y0^  reduce  to 


11  22  12  21  „ 
y = or.  , y = 0ro  , y = y =0, 


(9.17) 


where  for  a homogeneous  rod,  o^>On  are  cons'tan'ts'  The  Helmholtz  free  energy 
function  \Ji  can  be  expressed  in  the  form"'" 


* = tFi+  1»k  + *e+  *T  ’ 


(9.18) 


where  the  subscripts  F1,F2,E,T  attached  to  \|r  on  the  right-hand  side  of  (9.18) 
refer  to  the  four  modes  of  deformation  mentioned  above  and  we  recall  that  we 
adopt  the  values  (9*1 0)  for  D^,D  after  the  discussion  of  invariance.  The 
values  of  the  constituent  terms  in  (9-18)  are  listed  below  and  we  also  include 
the  notation  m.  ,m0,m  for  the  components  of  the  couple  m acting  over  any 
section  £=  constant  of  the  rod  where  ; 


m = m.  e. 
~ xid 


(9.19) 


All  response  functions  n^n^jP^ ,§,5a>k,k^  are  assumed  to  vanish  in  the 
reference  configuration.  The  identity  (9.15)  and  the  inequality  (9*16)  are 
also  used  in  obtaining  expressions  for 


In  (9.18),  as  well  as  the  rest  of  this  section,  we  use  the  same  symbol  for 
a function  and  its  value. 


Flexure  FI 


2p*Fl  k5V23 + k15K23 + k25K23e2 + k20^92J  ’ 

n23  = n2  k5V23  ’ ml  = P23  = k15K23  + ^k2502  ’ 

— “ ^2C^2  ” ^^25^23  y ?2  ~~  ci^2^  + c2®2  * *2  ” C2®2  9 ^2  ” ^2^2^^  * 


V23  = 623  + 5 K23  = ^23^  * 

^o/dC+Pfo  = Pd2uQ/Bt2  , 


(9.20) 


am^dC  - n2  + P^3  = ^2^ b 2 

pdi^/at  = p(s2+|2)  - ak£/a;  , 


Cg  ^ 0 , S 0 


Flexure  F2 


2p^F2  k6Yl3  + kl6K13  + k24Kl3  9l  + k19('0l'1  ’ 

nl3  = nl  = k6Yl3  5 m2  ~ ~ pi3  kl6K13  " ^k249l  » 

^19®1  ""  ^24*13  * ?i  = ^1^13  + ^2 91  ’ ?i  - ^29l  ’ ki  = ei^®i/^  > 


^13 + 

Bn^^/ac  + pfx  = pa2u1/at2  , 

Biig/a;  + n±-  pjt13  = - pQ'1d261o/at2  , 


(9.21) 


Pa\/at  = pu^)  - d^/BC  , 


b2  s 0 » ei  5 0 • 


1 


'T 


Extension  E 


2p*E  “ k1Y11  + k2Y22  k3Y33  + k7YHY22  + kBYllY33  + kQY22Y33 


+ k10Kll  + *11*22  + k17KHK22 


*21Y119  + k22Y229  + k23Y339  + kl8^0; 


nli  2klYll  + k7Y22  + k8Y33  + k219  ’ n22  = k7Yll  + 2k2Y22  + k9v33  + k229 


n3  k8Yll  + k9Y22  + 2k3Y33  + k23e  ’ 

P11  = k10Kll  + ^k17K22  ’ P22  = K7Kll  + kllK22  ’ 

5 = 0 > P'H  = - 2k2iVn  - 2k22Y22  ‘ ^k23Y33  " kl89  ’ k = sia0/aC 


Y11  2 ^11  5 Y22  2 ^22  ’ y33  , 


K-y_  - d&j^/aC  > ^22  = ®®22^^  ’ 


dn3/3C  + pf3  = pB2u3/at2  , 


2 2 

dPi;L/dC  + PAn  ' nn  = P®!3  6u/at  » 


2—  2 

ap22/3C  + P loo"^oo  = p aob  t , 


22  "22  K“2W  22' 
paVat  = pe  - ak/ac  > g!  ^ ° 


(9-22) 


^In  writing  the  coefficients  k-,  and  kg  an  overbar  is  used  to  avoid  confusion  with 
the  entropy  fluxes  k^,k2  in  FI  and  F2. 
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Torsion  T 


2ptT  +k12K12  + k13K21  + kl4IC12K21  ’ 

n12  = n21  = ^4^12^21^  ’ P12  ” P21  = m3  ’ 

p12  * k12K12+^kl4'C21  ’ P21  = k13K21  + 2kl4K12  ’ 

r^2  ~ ^21  ~ ^12  + ^21  * k12  = ^ 1 g^C  > = 36C)1/d^  > 


21  21' 


zp12/zc-n12  + pi12  = pa1a2612/at2  , 

2 2 

ap21/dC  - tt21  + pa21  = p<*2d  621/at 


(9.23) 


4o. 


11.  Straight  orthotropic  rod. 

V.'e  consider  the  special  case  in  which  the  direct  theory  models  the  main 
rr:perties  of  a straight  three-dimensional  linearly  elastic  rod  whose  line  of 
centroids  is  along  the  constant  unit  vector  e^.  Each  section  of  the  rod  is 
symmetric  with  respect  to  two  directions  specified  by  constant  unit  vectors 
e.  ,e  . In  addition  the  rod  is  orthotropic  with  respect  to  the  constant  ortho- 
r.crmal  basis  = (e^e^e^)  and  it  is  homogeneous  with  constant  density  p . 

Ihe  coordinates  along  the  basis  e.^  are  denoted  by  = (x,y,z),  the  constant 

— * 

reference  temperature  by  0 and  0 is  temperature  which  is  zero  in  the  reference 
* * 

state.  If  u =uiei  is  the  displacement  vector  and  e . the  linear  components  of 


the  strain,  then 


(10.1) 


where  ( ) . = d(  )/dx. . The  Helmholtz  free  energy  function  now  takes  the 
j1  1 

form 


P * 2 Ckieijek£  “ij'-ij 


cijeij9"-i  P*(c/e)0*2- p\e*  , 


(10.2) 


where 


ij  kjfc  ji  id 

c.  = c..  , c,  " = c. * = c.°.  = c_“  . 

id  di  k£  id  k£  jfck 


(10.3) 


The  coefficients  c >ck£  are  constants  and,  since  the  rod  is  orthotropic, 

the  only  nonzero  components  are+ 

11  11  11  22  22  33  23  13  12  . , 

cll  ’ c22  * c33  ’ Cll  ’ c22  ’ c33  ’ c22  ’ c33  » c33  ’ c23  ’ c13  * cl2  * (10-4) 

The  corresponding  components  of  the  stress  tensor  t and  the  entropy  T)  are 

^ J 

given  by 


+Eee  Green  and  Zerna  [9*  §5*1+]. 


Ul. 


\ 


(10.5) 


The  position  vector  of  points  in  the  rod  in  its  reference  state  is  given 
by 

R = x.e  = xe,  + ye  + ze  (10.9) 

~ 1~1  /-lL  r*CL 

1 2 

and  this  corresponds  to  (9«3)  with  x=£  , y=  £ , z=£  and  D1  = = 1.  The 

integrals  in  (6.6),  (6.7),  (6.8),  (6.9),  (7-8)  to  (7*13)  and  (7.16)  are  all 
with  respect  to  cross-sectional  areas  in  the  reference  configuration  of  the 
rod  and  their  boundaries 

F(x,y)  = 0 , (10.10) 

with  X.=  e^  and  \a  = 0 in  (6.10).  Then,  from  (6.6)  with  £ =x,  £ = y,  £-  = z,  we 
have 


* * 11  22  * 
n = p , X.  = p = Ap  , X»1  = \y  = p I2  , \a2  = \y  = p ^ , (10.11) 

! i 

L i 

* 

where 


I 


the  integrals  being  over  any  cross-section  of  the  rod. 

In  order  to  make  identification  of  the  thermoelastic  coefficients  occurring 
in  the  direct  linear  theory  in  (9.20)  to  (9.23)  we  make  some  use  of  the  linear 
results  (9-11)  and  (9-12),  together  with  the  approximations  (6.5)  and  (7.15). 
Thus , 


jl  = i+xIi  + y£> = K+^ii  + yWti  > 


e = e + x91  + ye2  , 


(10.13) 


Vij=6ij+V  ’ Si  = &ai/aC  ’ Kai  = a6oi/^  • 

We  consider  first  the  thermal  coefficients  in  (9.20)  to  (9.22)  whose  values  may 
be  identified  by  a direct  use  of  the  formulae  (7.8),  (7.11)  to  (7.13)  and  the 
results  (10.5)2,  (10.6),  (10.7)  and  (10. 13).  Thus 


:1  = ° » \ = c2  =-Ad22/(p0)  , e2=- I^/e  , 

*20  = ‘ p cI]/®  » k25  = ‘ 2I1C33  » 


(10.14) 


bL  = 0 , b2  = b = - Adn/(p0)  , eL  =-  I2d  /0 


k19=-pcI2/0  , k24=-2I2c33  , 


(10.15) 


*21  = - Aci;l  , *22  = ‘ Ac22  * k23  " Ac33  * kl8  = ‘ P Ac/e  , 


gl  = ' AdW 0 • 


(10.16) 


There  are  no  thermal  coefficients  in  (9.23)  for  the  torsion  group  T. 

Turning  to  the  mechanical  coefficients  k^...,^,  we  note  that  values  for 
many  of  these  have  been  given  for  isotropic  rods  by  Green  et  al.  [5],  Here 
we  follow  similar  procedures  but  now  applied  to  orthotropic  rods  which  also 


f 

A 


have  two  geometrical  axes  of  symmetry  in  each  cross-section.  We  consider  first 
the  flexure  groups  FI  and  F2  in  (9.20)  and  (9»2l).  A comparison  of  the  solutions 
of  equations  (9-20)  sr.d  (9-21)  for  the  problem  of  pure  bending  of  the  rod  by 
couples  over  its  ends,  with  the  corresponding  solution  of  the  three-dimensional 
equations  suggest  that  we  take 


: k15  Ii/S33  > ki6  I2//s33  ’ 


(10.17) 


where  s ^ are  the  inverse  coefficients  to  c1^  defined  by 


ij  mn  i/.m.n  , 0m,ns 

Vi'ij '5<8k6j+»A)  • . <10-18) 

where  6^  is  the  Kronecker  delta  and  s^  are  subject  to  the  same  symmetry 
J rs 

restrictions  as  c^  in  (10. 3).  When  the  rod  is  isotropic  the  values  (10. 17) 
reduce  to  those  given  by  Green  et  al.  (I967)  and  Green  et  al.  (I97^b).  In 
order  to  specify  values  for  the  remaining  coefficients  in  F1,F2,  we 

consider  the  complete  solution  of  the  static  isothermal  problem  in  F2  in  which 
the  rod  is  unloaded  along  its  major  surfaces.  Thus,  from  (9. 21),  when  f 1 = = 0, 

the  static  equations  may  be  integrated  in  the  form 


n1  = N , m2  = M - Nz  , 


(10.19) 


+ (-S-  - R)z+S 

Ks 

O 


where  N,M,R,S  are  constants.  This  represents  the  solution  of  the  flexure 
problem  in  which  a beam  of  length  l in  the  region  OSzSi  is  loaded  along  an 
axis  of  symmetry  in  the  x-direction  by  a load  N.  If  the  couple  is  zero  at 


z = j l and  the  rod  is  clamped  at  the  end  z = 0 so  that  u.  = 0,  6=0  there,  then 


that  in  order  to  find  a suitable  value  for  k,  we  must  consider  the  complete 
three-dimensional  flexure  problem  including  displacements.  The  flexure  of  a 
symmetric  orthotropic  rod  by  a force  N at  the  end  z = 1,  along  an  axis  of 
symmetry  in  the  x direction,  is  specified  by  displacements 


2I2  ' 33' 


/ 11  2 22  2 w \ 33  / i.  . __ 

U = - (s^x  - s y )(X-z)  + (— — ) + Pz  , 


Ms33 

^ 'lz- 


Ni 


22 


v*  = - (X-z)xy 

2 


w = - — {s~(lz  - \ z )x  + $(x,y)  + Lxy  } - Px 


Also, 


s23  £i  t „13  £l  , o 

23  SX2  13  ^ 


subject  to  the  boundary  condition 


(10.21) 


where  P is  a constant  representing  a rigid  body  rotation  and  L is  a constant 
given  by 

11  22 

-2L  = ^ + f^+4  . 

23  13  23 
23  13  23 


(10.22) 


(10.23) 


,j 

i !► 

ii 

l r 


i. 


«tf  * <l*i  s33>y2  ‘ i ° 


(10.24) 


over  the  surface  (10. 10),  where  (X,m)  are  direction  cosines  of  the  outward 
normal  to  this  surface.  Recalling  the  approximation  (10.13)  for  u , from  the 
comparison  of  (10.21),  _ with  (10.20)  we  again  obtain  the  identification  (10.17) 
for  k^g.  In  view  of  (10. 20)^  we  choose  the  rigid  rotation  P in  (10. 21)^,  and 
the  function  $,  which  is  even  in  y and  odd  in  x and  is  arbitrary  to  the  extent 
of  an  additive  constant,  so  that 


$(0,0)  = 0 , P = - N ■"  ^Vlg  • 


45. 


(10.25) 


Then,  comparison  of  (10. 20)^  and  (10.21)  at  the  end  z = 0 of  the  rod  for  small 


values  of  x,y  yields  P=N/kg  sc  that 


F2 


k6  " J2'’  dx 


(10.26; 


Analytical  solutions  for  the  flexure  problem  specified  by  (10.23)  and  the 
surface  condition  (10.24)  may  be  found  for  a number  of  simple  cross-sections 
of  the  rod.  Vie  merely  quote  the  final  value  obtained  for  kg  for  a circular 
cross-section  of  radius  R,  namely 


F2  : k,  = 


nR£<3i fiifS.) 


Similarly,  when  the  section  is  circular,  for  flexure  FI  we  have 


FI  : k,.  = 


5 ‘-gp-g  ♦*■!!♦■£>  ' 

When  the  rod  is  isotropic,  these  reduce  to  the  common  value 


(10.27) 


(10.26) 


k5  = kg  = 2^ttR2(1+v)/(3+2v)  , 


(10.29) 


where  v is  the  Poisson  ratio.  This  may  be  compared  with  a value 


6m,ttR2  (1+v  )2/  ( 7+l4v+8v2 ) 


found  by  Green  et  al.  [5]  by  a different  procedure  which  is  slightly  greater 
than  (10.29)  for  values  of  v in  the  range  0Sv5j. 

For  the  torsion  group  T we  compare  first  the  solution  of  equations  (9.23) 
with  the  Seint-Venant  torsion  problem  in  linear  three-dimensional  theory  in  the 
manner  described  by  Green  et  al.  [10].  In  the  present  context,  for  an  ortho- 
tropic rod,  this  leads  to  the  choice 


T : k12  k13  , 2k12  - ^ = fi  > 


(10.30) 


where  £ is  the  classical  torsional  rigidity  for  the  rod.  With  the  choice 


46. 


(IO.30),  the  isothermal  static  equations  in  (9-23 ) divide  themselves  into  two 
groups , namely 


m3  = (k]_2  " 2 ki4^^6i2  " > ^3/az  + p(^12-  l21)  = 0 (10.31) 


P12  + P21  = (k12  + 2 klU)9(  612  + 521)/&Z  ’ n12  = n21  = V612  + 621}  ’ 

[a(p12  + p21)/az]  - 2tt12  + p(a12  + a21)  = o . 


(10.32) 


We  now  consider  the  three-dimensional  solution  for  displacements  and  stresses 
given  by 


* * * 

= Lyz  , u2  = Lxz  , u3  = 0 , 

( 

*12  ■ 2oSLz  ’ *13  ■ ^ ■ *23  ■ • 

The  system  of  stresses  in  (10.33)  is  in  equilibrium  and  can  be  maintained  with 
suitable  surface  tractions.  Corresponding  to  this  we  see  from  (10.13)  that  an 
exact  solution  of  equations  (10. 31)  and  (10. 32)  is 


(10.33) 


®12  ®21  = Lz  ’ p^12  p^21  n12  n21 


P12  + P21  = 2(K12  + 2 klU)L  » "12  = n21  = 2k4Lz  ‘ 

With  the  help  of  (6.7),  (6.8)  and  (10. 33),  we  make  the  identifications 


(10. 3^) 


P12  + P21  ” (I1C13  + I2C23)L  ’ n12  2Ac12Lz  * 


(10.35) 


This  leads  us  to  choose 


1 ! k12  + i'kllt-i(Il0i3tV’il)  ' H * Ac12  • 


(10.36) 


A complete  set  of  torsional  coefficients  is  given  by  (IC.30)  and  (10. 36). 


47. 


Finally  we  consider  the  extension  problem  governed  by  equations  (9*22). 


Following  Green  et  al.  [5]>  the  values  cf  many  of  the  coefficients  can  be 
specified  by  a comparison  of  static  solutions  of  (9-22)  with  homogeneous 
deformations  of  the  three-dimensional  rod.  This  yields  the  values 


, 11 
Acll  » 

: , 22 

*2  ' “C22  ’ 

k3 

_ 1 
“ 4 

a 33 
Ac33 

A 11 

AC22  ’ 

, . , 11 
k5  = * Ac33  ’ 

k9 

_ 1 

Ac22 

33 

(10. 37  J 


When  the  rod  is  isotropic,  the  results  10. 37 j reduce  to  those  given  in  Green 
et  el.  [5].  Next,  we  consider  the  three-dimensional  solution  for  displace- 
ments and  stresses  given  by 


* 

U1  = 

* 

Lxz  , = Myz  , 

* 1 2 

u^  = J Nz 

tll 

, t = 0 
’ 12 

'22  = 

(C22L+C22M  + C33N)Z 

t = c13- 
* 13  C13 

'33  = 

(C33L+C33M  + C^N)Z 

t = c23l 
’ 23  23 

(10.38) 


where  L,M,N  are  constants.  The  system  of  stresses  in  (10. 38)  satisfies  the 
equations  of  equilibrium  if 


(c13  + c^)L+ (c^  + c^^M+c^N  = 0 , (10.39) 

and  can  be  maintained  by  suitable  surface  tractions.  Recalling  the  expressions 
(10.13),  it  follows  that  the  corresponding  results  in  the  direct  extensional 
theory  (9.22)  are: 


i 


V * 


Yu  = 26u  = 2Lz  , v22  = 2622  = 2Mz  , y33  = 23^/ dt  = 2Nz  , 

K11  = L ’ *22  = M ’ 

pAil  = nll  = a(ciil+C22M  + C33N^Z  5 

p£22  = n22  = a(C22L  + C22M  + C33N^Z  ’ 

pf3*A<c33Ltc33M+c|")  -°  ’ 

P11  = k10L+2  ki7M  ’ p22  = 5 k17M  + kllL  ' 

With  the  help  of  (6.7),  (6.8)  and  (IO.38)  we  make  the  identification 


(10.40; 


P11  I2C13L  ’ P22  " I1C23M  ' 


(10.41) 


Comparison  of  (10. 4l)  with  the  corresponding  expressions  in  (10. 40)  leads  us 


to  choose 


E : k10  " T2C13  ’ kll  = I1C23  * k17  ~ ° 


(10.42) 


for  the  values  of  k^,!  , ,k  . A complete  set  of  values  of  coefficients  for 
the  extension  problem  are  contained  in  (10.16),  (IO.37)  and  (10.42). 

When  the  rod  is  isotropic  and  has  circular  cross-sections  of  radius  R 
we  see  from  (10. 30),  (IO.36),  (10.37)  and  (10.42)  that 

kio  = kll  kl2  k13  = ¥ uTTR  * kl4  = k17  = 0 * 

k^  = 2k^  - k = ^nR2  . 


(10.43) 


Previously  (Green  et  al.  1974b),  all  constitutive  coefficients  in  the  case 
of  an  isotropic  rod  of  circular  cross-section  were  determined  except  k}_0,kll  and 
k^.  The  above  more  general  development,  which  is  valid  for  orthotropic  rods  of 
circular  cross-section,  upon  specialization  also  provides  appropriate  values  for 
the  coefficients  k^k^  and  k1?  in  the  isotropic  case. 


rT^^‘  y *+-  ■ 
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theory  or\traight,  elastic,  crthotropic  rods,  including  the  determination 
of  the  relevant  constitutive  coefficients. 
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